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Abstract 

We present high-precision non-relativistic variational calculations of bound 
vibrational-rotational state energies for the H2 + and D2 + molecular ions in 
each of the lowest electronic states of E s , E u , and H u symmetry. The calcula- 
tions are carried out including coupling between E and IT states but without 
using the Born-Oppenheimer or any adiabatic approximation. Convergence 
studies are presented which indicate that the resulting energies for low-lying 
levels are accurate to about 10~ 13 . Our procedure accounts naturally for the 
lambda-doubling of the il n state. 

PACS numbers: 31.15.Ar, 31.15.Pf, 33.15.Fm 
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I. INTRODUCTION 



There are many calculations of bound state energies of the hydrogen molecular ion H 2 + 
using the Born-Oppenheimer approximation or various adiabatic approximations and there 
are a number of studies that investigate deviations of energies from the Born-Oppenheimer 
values. The present work is a systematic high precision nonadiabatic[] study of H2 + and 
D 2 + in each of the lowest electronic states of E 9 , E u , and U u symmetry carried out using 
variational basis sets. It is motivated by recent precise experimental spectroscopy of Rydberg 
states of the hydrogen and deuterium molecules that has led to accurate experimental values 
of the the electric dipole polarizability of the corresponding molecular ions in their ground 
states [0J. These experiments were followed by several papers detailing various nonadiabatic 
calculations of the electric dipole polarizability |2|-|5|. The present paper is the first in a 
series. We are using the eigenstates studied in the present work in a study of the electric 
dipole sum rules for H 2 + and D 2 + , including the polarizability. 

Several investigators have performed nonadiabatic calculations on the ground electronic 
state of H 2 + since Hunter and Pritchard [|] and Kolos [|7| reported the first precision calcu- 
lations. The most accurate calculations used variational basis set methods MM, variation- 
perturbation methods [10 , 11 1 , and artificial channel scattering methods |T2|JT3[| . Variational 



basis set calculations can be in principle quite accurate but appear to have been applied 
only to the lowest-lying eigenvalues of the E 5 symmetry. The variation-perturbation and the 
artificial channel methods yield energies for all of the vibration-rotational levels and have 
been applied to the states of E 9 and E M symmetry. There are other approaches applied to the 
Eg symmetry that have not yet reported precision as great as those mentioned above such 



as the adaptive finite element method flU] , the generator coordinate method [15], quantum 
Monte Carlo [16|] and perturbative approaches [|17[]. Energy calculations up to 1980 were 
reviewed by Bishop and Cheung JIB] and a useful, more general review covering up to 1995 
can be found in ]19| . 



II. THEORY 



In this section we derive the Hamiltonian and introduce the basis sets we used. Other 
derivations can be found in Refs. [pO| - |23|JlO| i P^I ■ Some of the operators we use were introduced 
in those references and Ref. [25]. Our intention is to avoid writing explicit matrix elements 
until the last steps and the spirit of the present derivation is closest to the derivations in 
Refs. |2gj2|j2|]. 



A. Hamiltonian 

In a space-fixed frame and with the center of mass motion removed the Hamiltonian for 
the homonuclear one-electron diatomic molecule is 



1 We would prefer to use the term 'batic, which we coined to avoid the double negative implied in 
nonadiabatic, but clarity must yield to convention. 
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H = -iM-VI - [i + IM- 1 ] V 2 + V(v, R), (1) 

where 

and M = |M n , with M n the nuclear mass, r the position vector of the electron from the 
midpoint of the vector R joining the nuclei, and R = |R|. We use atomic units throughout. 
The electronic (cartesian) coordinates are to be held fixed in the space-fixed frame in carrying 
out the derivatives in the gradient operator V# appearing in Eq. (P f27| , |28fl . 



Following Ref. []29| we introduce the rotational angular momentum TZ implicitly express- 
ing the Hamiltonian in a rotating molecular fixed frame. The nuclear kinetic energy is 
written as 

2M ~ 2MR? [ dR K dR + K ) ' W 
Defining a rotational Hamiltonian 

tz 2 

H rot ~ (4) 

we write 

2M~ 2MR*dR dR + rot ' [) 

where the three spherical polar coordinates comprised of R and the two angles (contained in 
the TZ 2 operator of H TOt ) contain the information on the orientation of the molecular fixed 
frame with respect to the space fixed frame. 

Since here we are ignoring electron and nuclear spins, the total angular momentum is 
N = TZ + L, where L is the electronic angular momentum. Using TZ = N — L, we replace 
TZ 2 in Eq. @ giving 

H TOt = — ^(N - L) 2 = — ^(N 2 + L 2 - N~L + - N + L~ - 2N Z L Z ), (6) 
rot 2MR 2K ' 2MR 2K z zji 

where the superscripts on L + and L~ and subscript z on L z refer to the components in the 
molecule-fixed frame 



Changing the electron coordinates from cartesian to prolate spheroidal coordinates 
(A,/i, x), we have r = |r| = f(A 2 + /i 2 — l) 1 / 2 . The operator ^ in @ is taken with 
the electronic (prolate spheroidal) coordinates held fixed in the molecular fixed frame and 
can be expressed as 

d d , dr d , 
dR = dR > x ^~~dRlfr' ^ ' 

where the term ^ on the LHS of Eq. (|7]) refers to the derivative with the electronic (carte- 
sian) coordinates held fixed as in Eq. @. 
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Using the RHS of Eq. (^) in Eq. (|3|) we can write the kinetic energy operator as 

(8) 



V 2 

V R 



2M 2M 
where 

and 



d 2 2d 2Y d „ r 2 2 TT 



OR 2 RdR R 2 OR R 2 ' 



2 ^ ® 2 ® f Q s 

Y = r^- (10) 

and it is now understood that the electronic (prolate spheroidal) coordinates are held fixed 
where appropriate. 

We use the expression 

2 



^-^ = v 2 (11) 



r 



to combine Eq. (||) and (^|), yielding 

d 2 2d 27 a „ r 2 1 



2M 2M 



(12) 

We define for later use the coupling term 

1 



-N-L+ - N+L-) (13) 



2MR 2 

that enters from Eq. QT2D into the Hamiltonian. 

The potential energy is given in terms of the prolate spheroidal coordinate system 
0,/i,x) by 

K < A ^ B > = s-*(F^r (14) 

and the electronic kinetic energy operator by 

V 2 = (A/R 2 )[X + (A 2 - 1)^(1 - n 2 y l d 2 /d X % (15) 

where 

X = (A 2 - /i 2 )~ 1 [(9/9A)(A 2 - l)0/dA + (<9/<9/i)(l - /2 2 )d/dfj]. (16) 
and the operator Y, Eq. ([T0|) , becomes 

Y = (A 2 - ^ 2 ) -1 [A(A 2 - l)d/<9A + //(l - /i 2 )<9/c>4 (17) 



The terms in L can be reexpressed in the (A,/i, x) coordinates, see for example Ref. p0| . 

The remainder of the Hamiltonian derivation follows that of, for example M], and in 
this way the Hamiltonian reduces to effective matrix elements that may be evaluated as 
integrals over A, /i, and %. 
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B. Basis sets and trial functions 



For the electronic states of E g , S M and Tl u symmetry investigated here we used a basis 
set composed of functions of the form |24 



^ C P (A, fi, X) = (A 2 - 1) |A|/2 4 A '[«(A - 1)] exp[-| a (A - l)]Pi A| (/i) exp(*A X ), (18) 

with 6 = 0, ...,B and a a nonlinear parameter. We used values of A = —1,0, and 1. The 
values |A| =0 and 1 correspond, respectively, to £ and II states. For the £ 9 symmetry 
c = 0, 2, .., 2C and p = g, for the S u and U u symmetries c = 1, 3, 2C + 1 and p = u, and 
for the U g symmetry c = 2,4, 2C + 2 with p = g. 

The trial function for a particular set of states specified by A, p, and N has the form 

s 

R)= E ^ [M $^(A,/i,x)Xd(^) (19) 

s[6crf]=l 

where is given in Eq. (|1|) and where S = (B + 1)(C + 1)(D + 1). The index s = [bed] was 
filled in the order [{b, {c, {d}}}}, where {b}, for example, indicates a loop over all possible 
values of the index 6 = 0, B. The vibrational basis functions were of the form 

Xd {R) = (l/R^RY^L^R) exp(-i 7j R), (20) 

with d = 0, D. The vibrational state quantum numbers were identified with levels in the 
spectrum resulting from the diagonalization. The eigenvalues approach the exact eigenen- 
ergies behaving as expected by the Hylleraas-Undheim theorem [[H . 

Laguerre polynomials were used in the electronic basis because the integrals involved 
could be solved in closed form. Other possibilities explored such as Hermite polynomials did 
not offer this advantage. The electronic basis flT8]) is independent of R and is identical to that 



used by Moss and Sadler | 24fl . The vibrational basis is similar to theirs in functional form, 
but we used a different nonlinear parameter 7 that allowed us to avoid certain expressions 
involving hypergeometric series and thereby offered an apparent improvement in speed. We 
expect that the accuracy of our vibrational basis is at least equal to that of Moss and Sadler. 



III. CALCULATION 

Matrix elements of the Hamiltonian over the basis set functions and the overlap be- 
tween basis set functions were set up as four-dimensional integrals over A, n, x, and R. The 
evaluations reduce to integrals over A, /1, and R. The eigenvalues were obtained using the 
Rayleigh-Ritz method by solution of the generalized eigenvalue problem for the Hamiltonian 
and overlap matrices and iteratively varying the nonlinear parameters. Some details on the 
integrals and procedures are presented in this section. 



A. Evaluation of the integrals 

Consider the integrals over A and over R required for evaluation of the Hamiltonian and 
overlap matrix elements. Any integrals containing derivatives were manipulated to eliminate 
the derivatives by utilizing 
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^K(x) = -l£\(x) (21) 

and 

n 

K{x) = Y. L V\*) (22) 

fc=0 

to rewrite each integrand as a linear combination of integrals of the form 

/ dxx a+r L a m (x)L a n (x)e- x , (23) 

<J 

where r is an integer, r > 0. 

The resulting sets of integrals of form (^), and any other integrals of that form, were 
then manipulated to eliminate the powers of A. This was done by writing the product 
x r L'f n (x) as a linear combination of Laguerre polynomials with the same superscript. To 
this end, the expression 

xL a n (x) = (n + a)L a -\x) - (n + l)L a n ^(x), (24) 

derived using the summation definition for associated Laguerre polynomials, can be reduced 
using 

Lfa) = L?\x) - L£\{x) (25) 

to the desired expression, 

xL a n (x) = (2n + a + l)L a n (x) - (n + l)L a n+l (x) - (n + a^^x). (26) 



Substituting Eq. (|2q) into Eq. (|23|), each integral over A can now be expressed as a sum 
of integrals of the form 

/ dxx a L a m (x)L a n (x)e- x = 5 mn (m + a)\/m\. (27) 

The integrals involving fi could be performed through simple manipulations of associated 
Legendre polynomials. 

Coupling between states of different A introduced two problems. The first was that in 
order to carry out manipulations such as those used above leading to (|27|) , we required 
expressions for raising or lowering superscripts by more than unity. Using Eq. ( |22"D we 
derived the relation 

^) = E( i + * _1 )^to (28) 
and similarly from repeated application of Eq. ( p5]) we derived the relation 

K{x) = j:yi) k {^K + \{x). (29) 
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The second problem was the coupling of different 7 parameters. By using the same 
manipulations as for the A integral, we reduce the vibrational integral to a linear combination 
of functions /, where 



I(a, m, n, j it jj) = / dx x a L a m ('y { x)L*('y j x) exp(-|(7; + jj)x), 



(30) 



which can be reexpressed in terms of the hypergeometric function 2-^1 using Eq. (7.414.4) of 
Ref. IB2I as 



I(a, m, n, 7^ 7,-) = F(-m, -n; -m - n - a; 7 r 2 at ) 



[m + n + a)\ 



mini 



>a+l 



(-i) m 7-r m (7,+7,) 



where 



7rat = + 



-a-l 



(31) 



(32) 



The hypergeometric series terminates since m > and n > 0. Some additional notes on 
evaluating integrals of Laguerre and Legendre polynomials are given in [24]. Maple V was 
used to check the evaluation of the matrix elements and it was used to output them into 
Fortran code. 



B. Numerical procedures 

The trial functions (|T9|) have three sectors. They are comprised of two electronic sectors, 
labeled by the indices b and c and governed by the nonlinear parameter a, and one vibrational 
sector, labeled by the index d and governed by the nonlinear parameters (3 and 7. In our 
calculations each sector was treated separately in optimizing the nonlinear parameters and 
in studying convergence as the basis size was increased. The eigenvalues and wave functions 
were determined by solution of the secular equation using the lapack routines DSYGV 
and DSPGV, part of the math subroutine library dxml. The energy was further minimized 
by iteratively varying various nonlinear parameters (using a procedure discussed below) 
and rediagonalizing. For small basis set sizes we used a conjugate gradient method and 
then minimized by hand and for the larger basis set sizes we used an algorithm similar to 
Brent's [53|. Minimization of a was accomplished with standard algorithms. The optimum 



values for the parameters (3 and 7 were more difficult to determine for two reasons. First, 
(3 is integer and the necessarily discrete choices impeded the optimization; furthermore, a 
change in f3 does not correspond to a parabolic change in the value of the energy. Second, the 
nonlinear parameters j3 and 7 are intrinsically linked requiring simultaneous minimization. 

A general procedure was developed which allowed us to optimize a, /3, and 7 efficiently. 
Four steps can be identified. 1) We fixed (3 and 7 and then a was optimized for a minimum 
energy. 2) To minimize on j3 and 7 we fixed j3 and then minimized on 7. The parameter 
(3 was then varied by a large interval (about 6) and then we minimized again on 7. Some 
care was required in selecting what would be the optimum values of 7 as false local minima 
occasionally appeared. 3) Values of (3 within the final interval were searched for the optimum 
value with minimization on 7. 4) After all of the above a was reoptimized with the selected 
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/3 and 7. In all cases it was found in step 4) that the value of a was the same as that found 
in step 1), an important verification of our choice of final optimized nonlinear parameters. 

Having fixed the nonlinear parameters the basis set size was systematically increased to 
obtain precise eigenvalues by expanding each sector separately. Convergence to the final 
value was logarithmic. For H2 + in Figs. [I], 0, and [5] the convergence is demonstrated by 
plotting the difference between the energy for a particular basis set dimension and the 
energy for a basis set of dimension one unit larger. Results for D2 + are similar. For each 
figure, we begin with the final optimized wave function. The nonlinear parameters are not 
changed but the basis set dimension is set to B — 2, then index B is increased with the 
others held fixed at their optimized values and the difference between successive energies 
is plotted yielding the curves labeled U B (Electronic)" and similarly for C and D. For the 
S u states of H2 4 " and D 2 + convergence in the vibrational sector is slower than for the S 9 
and n„ states so we extrapolated to the desired numerical accuracy using linear regression 
on the log of the energy differences. Figure |2] illustrates the slow convergence but also the 
validity of the extrapolation. The basis set dimensions and nonlinear parameters for states 
with N = are given in Table | for S 9 symmetry in the first row under "Type I" and for 
S u symmetry in the first row under "Type II". 

For the states with N > 0, the off-diagonal term Eq. (|T^) in the Hamiltonian requires 
the inclusion of coupling between basis sets of S and II symmetry. Denoting the electronic 
basis sets by their value of A as |A) we set up matrix elements of the Hamiltonian using 
the rotated basis + | — 1)) and — | — !))• With it there is only coupling 

between |0) and ^(|+1) — |— 1)). A two by two matrix of matrices was created with the 
uncoupled Hamiltonian matrix elements for each basis set as the diagonal elements and the 
matrix elements of the coupling term Eq. (JT^) between the two basis sets as the off-diagonal 
elements. The energies of the states were determined by diagonalization of this matrix, 
while the energies corresponding to the uncoupled basis 4= (|+1) + | — 1)) were determined 
by diagonalization of the uncoupled Hamiltonian. For each state, the non-linear parameters 
and basis size were fixed at the values already determined for the minimum energies. Then 
the same technique used for the uncoupled energies was applied to the coupled basis sets to 
determine non-linear parameters and basis sizes that minimized the energy of the state under 
consideration. For example, when trying to determine the E u , v — 0, N — 1 energy, the S 
basis set parameters were held fixed at their uncoupled values, and the n basis set parameters 
were changed. The parameters for the coupling basis set were significantly different from 
those which minimized the energy in the uncoupled calculations, requiring six specialized 
parameters for each state when coupling was considered. The rate of convergence of the 
coupling terms is illustrated in Figs. |] and [5] for H 2 + . The energies converge logarithmically 
as each sector dimension is increased in turn. To evaluate the contribution of this small off- 
diagonal term to the energy many fewer basis set elements are needed than for the diagonal 
terms. The basis set dimensions and nonlinear parameters for states with N > are given 
in Table Q. For each symmetry there are two rows. The first row lists the dimensions and 
parameters for the primary symmetry used for all calculations and the second row lists the 
quantities for the additional symmetry required for N > entering through the coupling of 
Eq. (|TJ. 

The total number of basis functions used can be calculated from the data listed in Table | 
and is the sum of the values of S defined in Eq. (|l!|) entering for each symmetry. For example, 



S 



for H 2 + T, g , v = 0, N = 0, we used (13 + 1) x (5 + 1) x (13 + 1) = 1176 functions and for 
iV = 1 we used 1176 + (5 + 1) x (4+1) x (6 + 1) = 1386 functions. For H 2 + U u , v = 0, N = 1 
we used two runs, each corresponding to one of the rotated basis sets. For the uncoupled 
set, we had 910 functions, while for the coupled set, we used 910 + 270 = 1180 functions. 



IV. DISCUSSION 



Tables |H] an d [HJ compare the present calculations of nonadiabatic energies for H 2 + and 
D 2 + respectively with available precision calculations. In each table the vibration-rotation 
eigenvalues for the £ 5 symmetry are given first, followed by those for the £ u symmetry. 

For the £ 5 state the most precise variational basis set calculations are given for H 2 + in 
Refs. p4]- p6| . p!3| and for D 2 + in Refs. p^j35|jr3]| . Variation-perturbation calculations have 

and D 2 ~ 



for all the £ 9 vibration- 



been performed by Wolniewicz and Orlikowski |TT[ for H 2 ~ 
rotation states but the tabulated results include radiative and relativistic corrections and 
can not be compared directly with the present work. Using the artificial channel approach 
Moss carried out extensive nonadiabatic calculations of all the vibrational-rotational states 
,+ i r Jl cn^i r»„+ HJ2H f or the £ p states. His results with radiative and relativistic correc- 



of H 2 + f 131 and D 



tions are in good agreement with Wolniewicz and Orlikowski and he also presented energies 
without these corrections. In Tables [TT] and |TTT| the various calculations for the v = 0, N = 0, 
v = 0, N = 1, and v — 1, N — states are compared to our calculations. Results listed in 
Refs. fli~2|,|T3"f are converted from dissociation energies in wavenumbers to atomic units and 



combined with the asymptotic energy —M n /[2(1 + M n )]. Our results are consistent with 
and slightly improve upon the precision of previous calculations. 

Only a few high-precision calculations are available for the lowest states of £ u symmetry 
for H 2 + and D 2 + . Wolniewicz and Orlikowski used the variation-perturbation method and 
found 3 bound levels for H 2 + and 7 bound levels for D 2 + and gave energies of the levels with 
£-11 coupling included. Subsequently, Moss using the artificial channel method including 
£-11 coupling found results in agreement with those of Wolniewicz and Orlikowski for both 
H 2 + |~3| and D 2 + (T^j. Our £ u results are compared with these prior calculations in Tables [lj 
and |T|. For the v = 0, N = and v — 0, N = 1 states our energies are consistent with 
the others and of higher precision. However, for the D 2 + v — 1, N — state we found 
that a quite large basis set (B = 20, C = 11, D = 36 with a = 15.8, (3 = 37 and 7 = 2.6) 
was required to approach the energies given in Refs. [11,12]. Peek [37] showed that in the 
Born-Oppenheimer approximation the v — 1,N — vibrational wave function can have 
significant amplitude at values of R as large as several hundred a . Our electronic basis set 
is not explicitly dependent on R and this may account for the large basis size needed. Other 
methods |p~0V[T3|1 are based on coupled channel approaches that may be better at describing 
such diffuse vibrational states. 

There do not appear to be any published nonadiabatic energies for the lowest electronic 
state of n n symmetry of either H 2 + or D 2 + . Probably the most accurate study published 
is that of Bishop et al. 



38 



who investigated the U u energies of H 2 within the standard 
adiabatic approximation ||7|,|39|| . In Table [TV] the present nonadiabatic energies are compared 
to Born-Oppenheimer and standard adiabatic energies. The energy calculated in the Born- 
Oppenheimer approximation is a lower bound to the true energy while the standard adiabatic 
and nonadiabatic energies are upper bounds [5D,32[. The standard adiabatic energies were 
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calculated with the diagonal coupling of Ref. |38j rescaled to a proton mass of 1 836.152 701 
and the results differ in the seventh decimal place from the values reported in The 
present nonadiabatic results lie above the Born-Oppenheimer energy but below the standard 
adiabatic energy as expected [p2fl . 

The energies in Table [TV] were calculated without the consideration of Eq. (|I3"D leading to 
one level for each value of N. With the inclusion of the coupling term ( |13"D as described above 
in Sec. |111 B| our calculations exhibit lambda-doubling in the eigenvalues of II symmetry. In 
Table [V| calculated eigenvalues for the v = and 1 states with N — 1 are presented for H 2 + 
and D 2 + . For each value of v the first row gives the energy of the shifted level resulting from 
the diagonalization of the matrix coupling |0) and — |— 1)) and the second row gives 

the energy of the other, unshifted, level. The energy difference between the two levels is the 
lambda-doubling. 
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TABLES 



TABLE I. For rL + values of the dimensions B, C, and D and the optimized nonlinear pa- 
rameters a, (3, and 7. The values used for D2 + are identical except for the three values listed in 
parentheses. 









Dimension 




Nonlinear parameter 


Type 


Symmetry 


B 


C 


D 


a 


P 


7 


I 




13 


5 


13(17) 


3.1561 


67 


37.0 






5 


4 


6 


3.0 


79 


42.0 


II 




14(10) 


11(9) 


30 


15.8 


43 


3.1 




n u 


5 


5 


11 


13.0 


97 


7.4 


III 


n u 


9 


6 


12(19) 


6.0 


125 


16.5 






8 


5 


4 


5.0 


47 


3.86 
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TABLE II. Comparison of nonadiabatic vibration-rotation energies for H2 + for each of the 
lowest electronic states of £„ or Xl w symmetry. Calculations with TV > include the coupling term 
of Eq. (|lg|). Unless indicated otherwise all calculations correspond to a proton mass of 1 836.152 701 
in units of the electron mass. 



State 




Author (Year) 


Ref. 


Energy 


Eg,V = 0,N = 


: 


Bishop and Cheung (1977) a 


@ 


-0.597139 062 5 






Bishop and Solunac (1985) a 


H 


-0.597139 06318 






Moss (1993) 


H 

<\ _J 


-0.597139 0631 






Gremaud et al. (1998) 


@ 


-0.597139 063123(1) 






This work 




-0.597139 063123 9(5) 


Zg,V = 0,N = 


= 1 


Moss (1993) 


m 


-0.596 873 738 9 






This work 




-0.596 873 738 832 8(5) 


Vg,V=l,N = 


: 


Bishop and Cheung (1977) a 


m 


-0.587155 675 8 






Moss (1993) 


@ 


-0.587155 679 2 






flrPTnaiiH pt nl (1 QQ8^ 

v_Jl clllctLlvl CI dl. 


[M 


— D ^87 1 ^ R7Q 91 9(1 \ 






This work 




-0.587155 679 213 6(5) 


Z u ,v = 0,N = 


= 


Wolniewicz and Orlikowski (1991) 





-0.499 743 49 






Moss (1993) 


@ 


-0.499 743 502 2 






This work 




-0.499 743 502 21(1) 


Z u ,v = 0,N = 


= 1 


Wolniewicz and Orlikowski (1991) 


Q 


-0.499 739 25 






Moss (1993) 


@ 


-0.499 739 268 






This work 




-0.499 739 267 93(2) b 


a Proton mass 
b For this enerj 


1836.15 
5y, the b 


asis set had dimension B = 16. 
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TABLE III. Comparison of nonadiabatic vibration-rotation energies for D2 + for each of the 
lowest electronic states of £„ or X]^ symmetry. Calculations with iV > include the coupling 
term of Eq. (13). Unless indicated otherwise all calculations correspond to a deuteron mass of 
3 670.483 014 in units of the electron mass. 







State 




Author (Year) 


Ref. 


Energy 


y ' 


V 


= 0,N = 





Bishop and Cheung (1977) a 


Pt 

LL_JJ 


-0.598 788 782 










Bishop and Solunac (1985) a 


m 


-0.598 788 782 22 










Moss (1993) 


@ 


-0.598 788 784 










This work 




-0.598 788 784 330 8(1) 




V 


= 0,N = 


1 


Moss (1993) 


@ 


-0.598 654 8731 










This work 




-0.598 654 873 220 5(5) 




V 


= 1,N = 





Bishop and Cheung (1977) a 





-0.591603115 4 










Moss (1993) 


@ 


-0.591603122 










This work 




-0.591603121903 2(1) 




V 


= 0,N = 


o 


Wolniewicz and Orlikowski (1991) 


im 


-0.499 888 93 










Moss (1993) 


@ 


-0.499 888 937 5 










This work 




-0.499 888 937 71(1) 




V 


= 0,N = 


1 


Wolniewicz and Orlikowski (1991) 


Q 


-0.499 886 38 










Moss (1993) 


@ 


-0.499 886 382 5 










This work 




-0.499 886 382 63(1) 




V 


= 1,N = 





Wolniewicz and Orlikowski (1991) 





-0.499 865 21 










Moss (1993) 


M 


-0.499 865 2210 










This work 




-0.499 865 217(5) b 



a Deuteron mass 3670.48 

b For this energy, the basis set had dimensions B = 20, C = 11, D = 36 with nonlinear parameters 
a = 15.8, = 37, and 7 = 2.6 as discussed in the text. 
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TABLE IV. For H2 + the first several eigenvalues of the U u symmetry with N = 1 calculated 
nonadiabatically compared with Born-Oppenheimer and standard adiabatic calculations, respec- 
tively. For the present calculations, col. 4, the coupling term ( fLip has not been included. 



Vibrational state Born Oppenheimer Standard Adiabatic Present 3, 

-0.133 905 216 5 -0.133 841244 8 -0.133 841939 2 

1 -0.132 752 8516 -0.132 689153 4 -0.132 689 7691 

2 -0.131660 9817 -0.131597475 8 -0.131598133 6 

3 -0.130 6313519 -0.130 567 953 2 -0.130 568 676 9 

4 -0.129 666127 2 -0.129 602 748 3 -0.129 603 5416 



a Nonlinear parameters a = 6.0, (3 = 125, 7 = 16.5 with B = 9, C = 6, D = 24. 

TABLE V. Lambda-doubling in nonadiabatic vibration-rotation energies of H2 + and D 2 + for 
the lowest electronic state of H u symmetry for v = and 1, with N = 1. For each value of v the 



first row gives the energy of the shifted level arising from the coupling term in Eq. (13) and the 
second row gives the energy of the other, unshifted, level. 

Ion State Energy 

H2+ n„, i> = 0, iV = 1 -0.133 841940 395(5) 

-0.133 841939176 3(1) 
n u ,u = l,iV = l -0.132689769820(5) 

-0.132 689 7691218(1) 

D 2 + U u ,v = 0,N = 1 -0.134 052118 044(5) 

-0.134 052117 739 8(1) 
U u ,v = l,N = 1 -0.133 224 515 520(5) 

-0.133 224 515 448 7(1) 
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FIGURES 




"J ^ I i I i I i I i I i I i I 

2 4 6 8 10 12 14 

index n 

FIG. 1. Convergence study for the ground state S 9 energy of H2 + with v = 0, N = 0. The 
three basis sectors are fixed at their optimized dimensions for B, C, and D. Then for each sector, 
in turn, the index of the basis set B, C, or D, is set back to 2 and the value is increased until the 
optimized value of B, C, or D is reached again. Each line represents the log 10 of the energy for 
the index value n subtracted from the energy for the previous index value. (For sector B we have 
omitted the energy £12 •) 
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FIG. 2. Convergence study for the S u energy of H2 + with v = 0, N = 0. 
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index n 

FIG. 3. Convergence study for the Tl u energy of H2 + for the v = 0, N = 1 state with with no 
coupling to the S M symmetry included. 



17 



-5 



•- -• B (Electronic) 

♦ ♦ C (Electronic) 

■ ■ D (Vibrational) 



-10 




-15 




6 8 
index n 



10 



12 



FIG. 4. Convergence study for the energy of H2 + in the S u , v = 0, N = 1 state for the basis set 
of II M symmetry entering in the calculation. The S n symmetry basis set is fixed with the optimized 
size and nonlinear parameters listed in Table | for the calculations of this plot. 
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FIG. 5. Convergence study for the energy of H2 + in the Il n , v = 0, N = 1 state for the 
basis set of E u symmetry entering in the calculation. The Tl u symmetry basis set is fixed with the 
optimized size and nonlinear parameters listed in Table | for the calculations of this plot. 
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